Graphene on incommensurate substrates: 
trigonal warping and emerging Dirac cone replicas with halved group velocity 
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The adhesion of graphene on slightly lattice-mismatched surfaces, for instance of hexagonal boron 
nitride (hBN) or Ir(lll), gives rise to a complex landscape of sublattice symmetry-breaking poten- 
tials for the Dirac fermions. Whereas a gap at the Dirac point opens for perfectly lattice-matched 
graphene on hBN, we show that the small lattice incommensurability prevents the opening of this 
gap and rather leads to a renormalized Dirac dispersion with a trigonal warping. This warping 
breaks the effective time reversal symmetry in a single valley. On top of this a new set of massless 
Dirac fermions is generated, which are characterized by a group velocity that is half the one of 
pristine graphene. 
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Introduction - One of the main experimental challenges 
towards the realization of next-generation graphene elec- 
tronics technology is the possibility to access the low 
energy Dirac point physics. Silicon oxide (SiC^) sub- 
strates, for instance, are not ideal for graphene because 
of the trapped charges in the oxide. These impurity- 
induced charge traps limit the device performances and 
make the low energy physics inaccessible Q] . It has been 
recently shown that placing graphene on hexagonal boron 
nitride (hBN) yields improved device performances [2] - 
graphene on hBN can have mobilities and charge inho- 
mogeneities almost an order of magnitude better than 
graphene devices on SiC>2. 

hBN is interesting because it has the same honeycomb 
lattice structure of graphene, but only with two atoms in 
the unit cell, B and N, that are chemically inequivalent. 
Precisely this causes hBN to be a wide bandgap insula- 
tor. When graphene is placed on top of a hBN surface, 
the lowest energy stacking configuration has one set of C 
atoms on top of B and the other C sublattice in the mid- 
dle of the BN hexagons [31 H] - assuming perfect lattice 
matching between graphene and hBN. Consequently the 
substrate-induced potential breaks the graphene sublat- 
tice symmetry. This leads to a gap at the Dirac point 
and hence a robust mass for the Dirac fermions. First 
principles band structure calculations [3] put this gap at 
~ 50 meV - an energy roughly twice as large as fcgT 
at room temperature. However, recent scanning tunnel- 
ing microscopy experiments [5, 6 J do not detect a sizable 
bandgap. 

Within an effective continuum approach, here we show 
that this discrepancy originates from the 1.8 % lattice 
mismatch [7] between graphene and hBN which leads to 
a Moire superstructure with periodicity much larger than 
the graphene lattice constant. In this Moire lattice, car- 
bon atoms are embedded in a local environment of boron 
and nitrogen atoms that is varying continuously and pe- 
riodically. This leads to a complex landscape of local 
sublattice symmetry-breaking terms which prevent the 



opening of a band gap at the Dirac point. Due to the in- 
commensurability, the Dirac cones are instead preserved 
in renormalized form, with a threefold global symmetry 
due to a substrate-induced trigonal warping, which is in 
excellent agreement with experimental observations [8]. 
In addition we also show that a new set of massless Dirac 
fermions is generated at the corners of the supercell Bril- 
louin zone. These quasiparticles are characterized by a 
collinear group velocity vp which in the relevant weak 
coupling regime equals one half of v F , the Fermi veloc- 
ity in pristine graphene. As a set of these newly gener- 
ated massless Dirac fermions does not overlap in energy 
with any other states, gated or doped graphene triangu- 
lar Moire superlattices provide a clear way to probe these 
Dirac fermions. 

Before presenting the calculations that explicate these 
results, we wish to point out that very similar physics 
arises for graphene on incommensurate substrates other 
than hBN, in particular for the experimentally relevant 
Moire superlattices formed by graphene epitaxially grown 
on Ir(lll) surfaces [HE]- As the (111) surface iridium 
atoms form a triangular lattice, there are two distinct lo- 
cal Ir-C configurations with a high symmetry |10| . The 
first one occurs when a C atom is on top of Ir, situ- 
ating its three neighbors in throughs between Ir sites - 
the natural equivalence of the two triangular graphene 
sublattices is therefore broken. In the other high sym- 
metry Ir-C configuration, the honeycomb carbon ring is 
centered above an iridium atom. In this case the ef- 
fect of the Ir charges does not break the sublattice sym- 
metry and therefore, no gap should open at the Dirac 
point. While H decorated graphene/Ir(lll) superstruc- 
tures have been reported to give rise to absolute band 
gap openings HH] , recent angle-resolved photoemission 
spectroscopy data have rather shown an anisotropic be- 
havior of the massless Dirac fermions close to the Dirac 
points due to an enhanced trigonal warping [8 . 

A number of interesting theoretical predictions exist 
on graphene superlattices. It is known that external 
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one-dimensional periodic potentials can lead to a huge 
anisotropic renormalization of the electronic spectrum 
[TT| fT2] , emerging zero modes [13] and even to a Landau- 
like level spectrum as a result of the presence of extra 
Dirac points P3] . Dirac cone replicas at different k point 
in the Brillouin zone (BZ) have been reported also in 
triangular graphene superlattices [T5] as well as in bi- 
layer graphene superlattices [THl Hi]. These findings, 
however, rely on a description that disregards local sub- 
lattice symmetry-breaking terms, which are crucial when 
investigating the opening/closing of gaps in graphene on 
slightly incommensurate hBN or Ir(lll). In this Letter, 
we investigate the modification of the electronic spec- 
trum of graphene Moire superlattices taking explicitly 
into account these essential, slowly varying, sublattice 
symmetry-breaking terms in an effective continuum ap- 
proach. 

Effective-Hamiltonian - We start out by taking into 
account the interaction induced by the substrate charges 
as an external electrostatic potential for graphene's Dirac 
electrons. The potential has a triangular periodicity that 
coincides with the arrangements of the centers of the 
BN hexagons (Ir atoms at the 111 surface): V(r) = 
J2g ^G elG r where the G's are the reciprocal lattice vec- 
tors and Vg the corresponding amplitudes whose mag- 
nitudes depend on the modulus of G alone. In the 
following, we restrict the sum to the six wavevectors 
with equal magnitude G — 4n/(3as), G/G = (±1,0), 
(±C0S7r/3,±sin7r/3) , where \fias indicates the BN 
hexagons (Ir-Ir) distance. As the mismatch between as 
and the graphene carbon-carbon distance a is small, we 
can evaluate the effect of the substrate-induced electro- 
static potential on the two triangular graphene sublat- 
tices A/B as the sum of products of rapidly varying 

Jvl B A/B 

parts exp (iGsr -v-' ), r- 1 being the actual atomic 
positions and G$r rescaled wavevectors with magnitude 
Gsr = 47r/(3a), times slowly varying parts exp (zG • r) 
which we treat in the continuum limit. The G's are 
the rescaled "coarse-grained" wavevectors with magni- 
tude G = 47r|<5a|/3a 2 where 8a = as — a indicates the 
lattice mismatch which without loss of generality has 
been assumed positive. As a result, the effect of the 
substrate charges leads to an average external potential 
acting equally on each carbon atom and a mass term 
breaking the graphene sublattice symmetry given by 
V±{t) = [V a {t)±V b {v)}/2 = V /2Y,g [l±ex P z0 6 ] x 
exp(iG-r). The contribution of the rapidly varying 
parts of the potential are now encoded in the non-trivial 
phase factors cj)^ = G/G ■ S/a where S = — a(l,0) is the 
graphene nearest-neighbor vector. 

Since the large periodicity of the Moire super- 
structure prevents intervalley scattering, we can de- 
scribe the low-energy quasiparticles near the corners 
K± = |27r/(3a), ±27r/(3\/3a)} of the graphene hexag- 
onal Brillouin zone as 4-dimensional spinors "J" = 
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FIG. 1: (Color online) (a) Renormalization of the Fermi veloc- 
ity at the Dirac point as a function of Vb. The points are the 
results of the exact diagonalization of the low-energy Hamilto- 
nian whereas the continuous line is the analytical result from 
second-order perturbation theory, (b) Energy of the Dirac 
point as a function of the external potential amplitude Vo- 
(c),(d) Fermi lines for different values of the potential ampli- 
tude Vb close to the K+ (c) and the K_(d) valleys. Energies 
and wave vectors are in units of TivfG and G respectively. 

[Vfc + ,A)Vte + ,BjV'K-,B:V>K_,A]j characterizing the elec- 
tronic amplitudes on the two crystalline sublattices, with 
an effective Hamiltonian 

H = v° f t <g) p • a + V + (r)r ® <r Q + V_(r)r 2 ® a z . (1) 

Here we use direct products of Pauli matrices cr xyz , 
(T = 1 acting in the sublattice space and T xyz ,To = 
1 acting on the valley degree of freedom (K±). For 
V±(r) = 0, H has a chiral symmetry which can be ex- 
pressed as |"H,To ® cr z \ = 0. This anticommutation 

relation implies that in each valley any eigenstate \l/ e 
with energy e has a particle-hole partner tq <S> a z ^ with 
energy — e. This property implies the doubly degener- 
acy of the zero energy states in each valley. In the 
presence of substrate-induced interactions of the form 
as in Eq. [T] the system still possesses a chiral symme- 
try provided the external superlattice potentials satisfy 
V±(r + T) = — V±(r). In this case it is possible to de- 
fine a new chiral operator [14] To <8> a z S where S is a 
shift operator S^r) = <I/(r + T). For the electrostatic 
potentials defined above, a translation vector for which 
the triangular potential V+(r + T) = — V+(r) is absent 
thereby implying particle-hole symmetry breaking and a 
consequent lifting of the zero energy states degeneracy. 



This, however, does not lead to the opening of any ab- 
solute band gap since the G = component of the local 
sublattice symmetry breaking terms identically vanishes. 
Therefore, the Dirac cones are preserved with the effect 
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FIG. 2: (Color online) Energy dispersion relations for a graphene triangual Moire superlattice with Vo = O.lh v%G close to 
the two graphene valleys K + (a) and K (b). The arrows indicate the corners K± of the SBZ where the new Dirac cones are 
generated . (c) Topology of the Fermi lines close to the emergent Dirac cones. Units are the same as in Figll] 



of particle-hole asymmetry eventually leading to a shift 
of the conical points [shown in FigQb)] of the two val- 
leys reminiscent of the graphene doping caused by ad- 
sorption of metal substrates [18] . We also find the Dirac 
cones to be renormalized in triangular Moire superlat- 
tices. In Figflja) we show the behavior of the collinear 
Fermi velocity at the conical points for different values of 
the interaction strength Vq . The substrate- induced inter- 
action leads to a decrease of the Fermi velocity as can be 
found in the weak potential limit by treating the effect 
of the electrostatic potentials in second-order perturba- 
tion theory [c.f. continuous line in FigfTJa)] according to 
which 



v F 



h 2 v° F 2 G 2 



The local sublattice symmetry-breaking term breaks 
the effective time-reversal symmetry on a single valley 
[15] , T = i{tq ® (j y )C 1 with C the operator of complex 
conjugation and T [ipK + ,A, ipK+,B, ^k_,b> V>k_,a] = 
[Vk + ,bi-^k±,A>^js.-,A>-$K-,b\- This is clearly vis- 
ible in Figfijc),(d) where we show the topology of 
the Fermi lines close to the Dirac points in the two 
graphene valleys. There is a trigonal warping which 
breaks the k — > — k symmetry of the Fermi lines, i.e. 
e(K±,k) 7^ e(K±,— k), consistent with the threefold 
symmetry of the bandstructure experimentally detected 
in Ir(lll) superlattices [8]. The trigonal warping has 
an opposite effect on the two valleys since the exter- 
nal electrostatic potentials do not break the true time- 
reversal symmetry interchanging the valleys [Tp T = 



(T y (g>a y )C with T[V'k + ,A,'0K + ,b,V'K_,b,^K_ 



lines fulfill e(K±,k) = e(K T , — k) as can been shown in 
the weak potential limit where the energy dispersion of 
the low-energy quasiparticles reads 



"V'k B'V'k a\- Hence, the Fermi 



e(K±,k) ~ hv F s\k\ 



3s |fc| 2 



1 ± V 3 sin 36*i 



(2) 



Here 9k is the angle of the vector k with respect to the k x 
direction and s = ±1 for quasi-electrons and quasi-holes 
respectively. By increasing the strength of the poten- 
tial Vq [c.f. Figfijc),(d)], the warping effect is enhanced 
and results in an anisotropy much larger than the one 
expected in freestanding graphene |S]. 

Emergence of Dirac cone replicas - We have also ob- 
tained the energy dispersion in the full supercell Brillouin 
zone (SBZ) by exact diagonalization of the Hamiltonian 
Eq. [TJ The effect of the substrate-induced external po- 
tential V± (r) has been incorporated into our calculations 
through the scattering matrix elements between the chi- 
ral eigenstates of the graphene quasiparticles 



<MK±,k) 



1pK±,A 
"0K±,B 



±( 



1 

±< ( 



A k-r 



Fig(2] shows the ensuing energy dispersion of the first and 
second bands above and below the original Dirac points 
in each valley. Contrary to Ref. [121 we do not find a gen- 
eration of Dirac cones at the six M points of the SBZ. 
Indeed, the energy separation between the first and the 
second band above and below the original Dirac points 
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goes to zero respectively at the (K±, K±), (K±, K T ) cor- 
ners of the SBZ. This qualitative difference is caused by 
the sublattice symmetry-breaking term V_ (r) which does 
not allow for a sixfold symmetry of the bandstructure. 
The topology of the Fermi lines close to the SBZ corners 
clearly shows the emergence of Dirac cone replicas. How- 
ever, while above the original Dirac points [c.f. Fig{2]jc)] 
these new massless quasiparticles are obscured by other 
states, below the original Dirac points there is an energy 
window where there are no other states than the new 
massless Dirac fermions. Therefore there is one energy 
value - apart from the original Dirac point - where the 
density of states (DOS) vanishes linearly. It is worth 
noticing that this asymmetry of the DOS reflects the 
particle-hole symmetry breaking discussed above. 

Further insight into the properties of the Dirac cone 
replicas is obtained by introducing an effective Hamil- 
tonian close to the three equivalent corners of the SBZ 
[20] , In the following we will restrict to consider the 
behavior close to the (K±,K T ) points, relevant for the 
Dirac cone replicas generated below the original Dirac 
points. In the absence of external electrostatic poten- 
tials V±(r) = 0, there are three degenerate hole excita- 
tions with energy e(K±,K T ) = —hvpK. This degener- 
acy is lifted by the substrate-induced electrostatic poten- 
tials and, as a result, one finds a singlet excitation with 
energy es(K±, K T ) = —hvpK — Vq and a doubly degen- 
erate state with energy e£>(K±,K T ) = —hvpK + Vq/2. 
It can be easily shown that the effective Hamiltonian in 
the vicinity of this doubly degenerate state corresponds 
precisely to a massless two-dimensional Dirac equation 
with Fermi velocity Vp = v F /2 and an isotropic disper- 
sion relation 



e D (K ± ,Sk ± ) = e D (K ±1 K T ) 



(3) 



where we introduced 5\i± = k — Kzp and s' = ±1 is the 
band index. The foregoing weak coupling analysis is in 
perfect agreement with the numerical results obtained by 
exact diagonalization of the Hamiltonian. This is shown 
in Fig(3^a),(b) where we plot the behavior of the Dirac 
point energy ei)(K±,Kzp) and the Fermi velocity v F as 
a function of the potential strength Vq. Away from but 
close to the Dirac points, we find a trigonal warping re- 
specting the threefold symmetry of the band structure as 
it is shown in Figj3]^c),(d) where we plot the Fermi lines 
for different values of Vq. 

Conclusions - By setting up an effective continuum 
approach, we have derived the electronic properties 
of graphene Moire superlattices generated by adhesion 
of graphene sheets onto lattice mismatched substrates. 
While the complex landscape of sublattice symmetry- 
breaking terms prevents the opening of a bandgap at 
the Dirac point, we have demonstrated that with mis- 
matched substrates one can tailor the low-energy band 




FIG. 3: (Color online) Behavior of the Dirac point energy 
en(K±,K T ) (a) as a function of the external potential am- 
plitude Vo- The continuous line is the result of the degenerate 
perturbation theory whereas the points are the result of the 
exact diagonalization. (b) Group velocity at the Dirac point 
replicas as function of the strength of the potential Vo- (c),(d) 
Fermi lines close to the emergent Dirac points e_o(K+,K_) 
(c) and £.d(K_, K+) for different values of the amplitude Vq. 



dispersion. In agreement with recent experiments [5] , we 
have found a threefold-symmetry of the band structure 
associated to a substrate-induced trigonal warping of the 
Dirac cones. In addition a new set of Dirac fermions is 
generated in graphene Moire superlattices. By properly 
and explicitly accounting for local sub-lattice symmetry 
terms, we have shown that these new quasiparticles are 
generated at the corners of the supercell Brillouin zone 
and are characterized by a collinear group velocity at the 
conical points of ~ v F /2. 
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